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Abstract 

Adiabatic quantum computation provides an alternative approach 
to quantum computation using a time-dependent Hamiltonian. The 
time evolution of entanglement during the adiabatic quantum search 
algorithm is studied, and its relevance as a resource is discussed. 
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1 Introduction 

How does quantum computation work? What are the necessary physical 
resources? And what can we learn from it about fundamental physics? 
Adiabatic quantum computation may provide some contributions to these 
qüestions. 

Quantum information science has gone through an amazing develop- 
ment in the past decade, and fascinating applications have been developed, 
including quantum teleportation or quantum cryptography. It is hoped that 
quantum computing will provide improved performance as compared to a 
classical computer, that it will provide an alternative route where miniatur- 
ization reaches a limit, and that it might even allow to solve problems that 



'based on a talk by DA given at Theory Canadà I, UBC, Vancouver, June 2005 



1 



a classical computer cannot solve. But despite the amazing success of quan- 
tum information science with regard to applications, there is still no general 
agreement about the fundamental physics principies or physical resources 
that make quantum computing (at least theoretically) work and that pro- 
vide a speed-up over classical computation. Candidates are entanglement 
or nonlocality, or the Hilbert space structure including the superposition 
principle, which suggests that some mechanism of parallel computing might 
be at work. 

In this paper, I study the question of resources within an alternative 
model of quantum computation, so called adiabatic quantum computation 
(see, e.g., 0). The equivalence of adiabatic and Standard model quantum 
computation has recently been shown in |2J. The advantage for the investi- 
gation of entanglement is that in the adiabatic approach, the computation is 
done through a time-dependent Hamiltonian, which allows to study the evo- 
lution of physical quantities like the entropy of entanglement as a function 
of time. Below, I will provide a brief introduction to entanglement and to 
adiabatic quantum computation and its realization of the search algorithm. 

1 will describe a method for calculating entanglement during the running 
time of the algorithm and present results for the original algorithm as well 
as modifications. Note that the presentation is aimed at a physics audience 
and does not require previous knowledge of quantum computation. 

2 Physical Resources in quantum information sci- 
ence 

The subfield of quantum information science that has evolved most rapidly 
and already seems to be moving from purely scientific research to indus- 
trial applications (e.g., 4 ) is quantum cryptography. In this case, as in the 
related quantum teleportation (see, e.g., jü] for an illuminating discussion), 
entanglement is used as a resource, and the measurement process is part of 
the protocol. 

In comparison to these areas, the progress in quantum computation has 
been more modest. Algorithms have been found for the testing of functions 
(Deutsch-Josza's algorithm jüj), for factorizing large integers in polynomial 
time (Shor's algorithm jjj), and for searching (Grover's algorithm 8 ). Al- 
though the efficiency of these algorithms is well established, the question of 
which physical resources allow them to work is still a controversial issue (see 
e.g. [S] for a discussion about possible implementations of these algorithms 
using Nuclear Magnètic Resonance and the question whether entanglement 
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has been involved). It is generally believed that entanglement is necessary 
for quantum computing as a nonclassical feature providing the computa- 
tional speedup ^ü]. It has also been suggested to consider this resource 
not as an intrinsic property of quantum physics itself, but something that 
is tied to a particular choice of mathematical formalism |12| . On the other 
hand, examples have been presented where quantum computing without en- 
tanglement is better than the classical version jllj . The search algorithm 
has even been implemented using the wave nature of classical Fourier òp- 
tics ^31- Summarizing, it seems that entanglement somehow plays a role in 
quantum algorithms. But it is not understood yet whether the entanglement 
is necessary for the speed-up, or whether it is created by the algorithm as a 
by-product, if this depends on the algorithm considered or if it is a general 
property of quantum computing, and how this is related to the Hilbert space 
structure or the formalism in general and other physical resources involved, 
like energy. 

2.1 Entanglement 

Entanglement as a non-classical physical phenomenon (non-classical corre- 
lations) has been known since the early days of quantum mechanics. It plays 
the leading role in the EPR paradox, although in that context, it is better 
known as non-locality - and was nick-named "spooky actions at a distance". 
When Bell established his famous inequalities, the discussion was moved 
from more philosophical considerations to experimentally testable physics. 
And yet, after all this time, a precise definition and especially a general 
measure for entanglement is still missing. Consider, for example, one of the 
Bell states, 1/V2(|00) + |11)). It cannot be written as a pro duet of single 
qubit states, in other words, it is not separable. Unfortunately, this version 
of stating entanglement is not very practical in applications. It is desirable 
to have a measure of entanglement that is easy to calculate, and that allows 
to compare different states in their degree of entanglement. For bipartite 
systems (i.e. systems consisting of two subsystems, A and B) in pure states, 
the von Neumann entropy of either of the subsystems is generally accepted 
as the measure of entanglement. It gives the amount of information of one 
qubit that can be obtained by making a measurement on the other qubit of 
a pair, or, in other words, it measures how complete the information about 
the subsystem is, using the density operators of the subsystems: 

E(pab) = S(p A ) = -Tr(p A \og p A ) = S(p B ) = -Tr(p B \ogp B ) (1) 
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with the traces over either subsystem pa = Ttb^Pab) and ps = TtaÍPab)- 
For a system in an entangled state, the state of the whole system is com- 
pletely known, but the state of each subsystem is not, i.e. it is in a mixed 
state. For a maximally entangled state of two qubits, e.g. a Bell state, 
Pa = 1/2/, and E(pab) = 1- K the subsystem is in a pure state, the en- 
tropy of entanglement for the composite system vanishes. The von Neumann 
entropy can be calculated from the (non-zero) eigenvalues A n of the reduced 
density matrix as 

5 = -Tr(plogp) = -^A n logA ri . (2) 

n 

For more general systems, for example for mixed states, several other, in- 
equivalent measures of entanglement have been suggested, but no generally 
accepted one has been found yet (see, e.g., ^3] for a discussion). 

3 Adiabatic Quantum Computation 

The building blocks of a quantum computer are qubits, quantum mechanical 
systems with a two-dimensional Hilbert space, which can be realized as, e.g., 
spin 1/2 partides. (Note that the general state of a qubit is a superposition 
of the two basis states, |<3?) = a|0) + (3\1), whereas a classical bit is either 
in state or state 1.) In the traditional model of quantum computation, 
quantum circuits are considered which consist of a series of unitary operators 
(gates) acting on single or múltiple qubits. At the end of the computation, 
a measurement provides the result. 

The idea behind adiabatic quantum computation is to use a more "natu- 
ral" description for the time evolution of the system, a Schròdinger equation. 
The computation is understood as a continuous time evolution from an ini- 
tial to a desired final state: The initial state l^o) is the ground state of a 
Hamiltonian Hq and (assumed to be) easy to build, whereas the final state 
l^i) is the ground state of a Hamiltonian H\ and encodes the solution of 
the computational problem. The time dependent Hamiltonian drives the 
system from the initial to the final instantaneous ground state. It is usually 
constructed clS cl linear combination 

H(t) = f(t)H + g(t)H 1 (3) 

with the monotonic functions /, g fulfilling /(0) = 1,/(T) = and g(0) = 
0,g(T) = 1. The system is supposed to stay in its instantaneous ground 
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state during this time evolution. The condition for this behaviour is pro- 
vided by the adiabatic theorem: Consider a system evolving according to a 
Schròdinger equation with a time dependent Hamiltonian. Its instantaneous 
energy eigenstates obey 

H(t)\E k ;t) = E k (t)\E k ;t). (4) 

If the system is prepared in its ground state l-E^O) at t = 0, after time 
T the probability for staying in the instantaneous ground state (up to a 
phase) is \(Eq;T\^(T))\ 2 > 1 — e 2 , provided that the Hamiltonian varies 
slowly enough, 



max <í<T 



Ei)t 



du 



dl 
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with the minimum gap between the two lowest eigenvalues 

9min = mm^i (í) - E Q (t)] . (6) 

Accordingly, the running time of the algorithm is given by the minimum 
time T m i n for the evolution of the system to be adiabatic, eq. (JïïJ). After this 
time, a measurement will give the desired solution with probability P ~ 1, 
i.e. with accuracy e given by eq.©. 



4 Adiabatic Quantum Search 

The search problem consists of finding a marked object (the "needle") in 
an unsorted database (the "haystack") of N elements with the minimum 
amount of computational work. A classical search would simply examine the 
items one by one and find the marked one after O(N) steps. Grover's original 
algorithm 8 allows to search an unstructured database quadratically faster 
than any classical algorithm. In the quantum version, the items of the 
database are represented by the computational basis states, and the haystack 
is given by an equally weighted superposition of these states, 

1 N ^ 

!*.) = ^ g 10. (T) 

i.e. each one has the same probability. Grover's algorithm works by in- 
creasing the amplitude of the marked element while decreasing all other 
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amplitudes, i.e. by rotating the state of the system towards the marked 
state. Its action on a state can be written as 

G|*> = (2|*)(*| - /)(/ - 2\m)(m\)\V) = (1 - -^)|*> + ^=|m>. (8) 



Note that the factor l/y/N is the overlap of the state \^f) with the marked 
state |m). After ko ~ \\ÍN iterations, the probability for obtaining the 
marked state after a measurement is P > 1 — 1/N, i.e. the complexity of 
Grover's algorithm is 0(y/N), a quadratic improvement over the classical 
algorithm. 

For the adiabatic quantum search, a Hamiltonian has to be constructed 
which has the corresponding effect: H(t) = f(t)Hç, + g(t)H\ with the initial 
Hamiltonian 

flo = I-|*o><*o|, (9) 
which has the "haystack" as its ground state, and the final Hamiltonian 

H 1 = I -\m){m\, (10) 

which has the marked state \m) as its ground state. In its simplest version, 
with linear functions f(t) = 1 — t/T and g(t) = t/T, the adiabatic search 
has a running time of O(N), i.e. no speed-up over the classical search pQ. 
This result can be improved by applying the adiabaticity condition locally, 
with respect to time, instead of globally ( to the entire time interval): The 
Hamiltonian can change faster when the gap between ground state and first 
excited state is large, which is the case at the beginning and the end of 
the evolution. Only around g m %n h needs to change slowly. Adjusting the 
evolution rate with nonlinear functions /(£) and g(t) gives a running time 
of T m i n = 0{^/~N)/e [Ej) corresponding to Grover's result. Recently it has 
been studied how a temporary increase in energy for the oracle term, given 
by the function g(t), can increase the speed-up further. For a large number 
n of qubits, the time T for the computation is bounded below by ^ü] 



ni <>m>— dl) 

where k = 0(1). This result suggests that the action associated with the 
oracle term can be considered a resource required for adiabatic search. In 
the next section, I will study if entanglement can be considered a resource 
for adiabatic quantum search as well. 
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5 Entanglement in a two-qubit toy model 



For the investigation of entanglement in adiabatic quantum search, I con- 
sider the simplest "quantum computer" possible, a system of two qubits. 
The computational basis states for this composite system are written as 
|0) ® |0) = |0)|0) = |00), |01), |10), |11). The initial state is 

|$(0)) = -^(|00) + |01) + |10) + |11)), (12) 

and the marked state is taken to be 1 00) . From these states, the time- 
dependent Hamiltonian is constructed as described in the previous section. 
Since for a quantum computer of this size there is no speed-up over classical 
computation, I can also use the simplest time dependence of the functions 
f(t) and g{t). From the Hamiltonian, the lowest energy state is calculated 
to be 

|£L(í)> = 1 = U4^L " 1), 1, 1, , (13) 



where £L = 1/2(1 - y/3(t/T) 2 - 3t/T+ 1) is the lowest energy eigenvalue. 
Note that this corresponds to the general result (9) in |16j for the case of 
n = 2 qubits and the functions f(t) = 1 — t/T and g{t) — t/T. In the 
adiabatic approximation considered here, the systems stays in this instanta- 
neous ground state throughout its evolution, which allows one to calculate 
the entropy of entanglement in a straightforward way: from the coordinates 
of the state in the computational basis, co(í), c\(t), C2(í), c^{t), I calculate 
the eigenvalues of its reduced density matrix as 

1 



M± (í) = -(l±Vl-4C2(í)) (14) 

with C 2 {t) = |co(í)c3(í) — Ci(í)c2(í)| 2 , the concurrence as a function of time. 
The entropy of entanglement for the ground state as a function of time is 
calculated from these eigenvalues as £ = — (/x + log/i + + /<í_log//_). The 
result is shown in Fig. ^(left)- Although both the initial and the final state 
are separable, the entanglement is non-zero throughout the running time of 
the algorithm. Is this entanglement necessary for the algorithm to work, and 
to provide the speedup, or is it just a byproduct of the time evolution? For 
a more systematic approach, I vary the initial state, thereby modifying the 
search to one in which not all items are present with the same probability. 
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Figure 1: Entropy of entanglement for the equally weighted superposition 
as initial state (left) and for a Bell state as initial state (right). 



E(t) E(t) 




Figure 2: Entropy of entanglement for two non-maximally entangled states 
as initial states. 

(Note that the overlap of the initial and the marked state has to be nonzero 
for the search to succeed). A different initial state has, in general, a value for 
the entanglement different from the equally weighted superposition. It also 
leads to a different Hamiltonian, which induces different interactions in the 
system and therefore a different entanglement during the algorithm. Fig. 
121 shows two typical results for the entropy of entanglement for two initial 
states with non-maximal entanglement. In either case, the entanglement 
reaches its maximum value of one during the evolution before going down to 
zero. In the second case, the initial state yields a Hamiltonian with the same 
spectrum as the equally weighted superposition, and the same Trnin\ only the 
initial entanglement is different. This would support the assumption that 
the running time of the algorithm is correlated with the spectrum of the 
Hamiltonian, as the adiabatic theorem suggests, and not necessarily with 
the entanglement. The entanglement for a search starting in a Bell state is 
shown in Fig. ^ (right), starting at the maximum value one and going down 
to zero for the final state. 



These four cases represent the results found when varying the initial 
state systematically. The entanglement either approaches the value zero 
monotonically, which seems to be the case for initial states with cq > C3, 
i.e. states "close" to the marked state, or it reaches a maximum during 
the running time. The value of this maximum seems to be independent of 
the initial entanglement, but the difference between maximum and initial 
entanglement seems to be larger for states with C3 3> cq. This would imply 
that for states "further away" from the marked state, a larger amount of 
entanglement is created during the search. 

In the cases studied, no correlation between the minimum energy gap 
g m in and the initial entanglement could be found. Since the gap is corre- 
lated with the running time through the adiabatic theorem, this would also 
support the conclusion that entanglement is not the crucial resource for this 
algorithm. 

6 Discussion 

From the results presented in the previous section, no clear statement about 
the relevance of entanglement for the adiabatic quantum search algorithm 
can be made. This is not too surprising, since we know that the speed-up 
in Grover's search algorithm comes from the overlap of the initial and the 
marked state, 1/N = (ip\m). The relevance of this overlap, or the fidelity, 
for the algorithm should be investigated further [2Ü]. Since iV is given by 
the dimension of the Hilbert space, it seems plausible that the Hilbert space 
structure is a crucial resource for the search algorithm, allowing a kind of 
paral·lel computation. But a caveat seems necessary here: The Hilbert space 
dimension is related to the dimension of the maximally entangled subspace 
[T7| . and so, in a way, to the number of available entangled states. This could 
explain why there seems to be a contradiction between the results obtained 
here for a small system, and the results for large N found in ^H] for the 
Standard Grover search: the more entangled the initial state is, the less well 
the algorithm performs. In that paper, the quality of performance is not 
the running time, but the probability for actually finding the marked state, 
which might not be comparable with the adiabatic running time. But the 
study of entanglement for systems with large N in the adiabatic approach 
is certainly a future goal and will partially be addressed in [SO]- We plan 
to report on a larger, more systematic study of entanglement and fidelity 
during the search and other algorithms in j^ü]- Using numerical methods 
to obtain more extensive data, we will also attempt a more detailed inter- 
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pretation of the resources used in adiabatic quantum computation. From 
studying structured adiabatic quantum search ^ü] for large N, we know that 
a more structured, i.e. more local Hamiltonian has a shorter running time 
than a Hamiltonian including interactions between all subsystems. Since 
entanglement is created (and destroyed) by interactions, this would support 
the assumption that for large systems, entanglement slows down the algo- 
rithm. A better understanding of entanglement would generally be helpful 
for the understanding of the dynamical evolution of composite, interacting 
many-particle systems. 
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